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Abstract 

We show that using the concept of the two-dimensional level number N± one can experimentally 
study of the nodal domains in a three-dimensional (3D) microwave chaotic rough billiard with the 
translational symmetry. Nodal domains are regions where a wave function has a definite sign. We 
found the dependence of the number of nodal domains H n ± lying on the cross-sectional planes of 
the cavity on the two-dimensional level number iVj_. We demonstrate that in the limit iV|_ — > oo 
the least squares fit of the experimental data reveals the asymptotic ratio Njv x /-^Y± ~ 0.059 ±0.029 
that is close to the theoretical prediction Hn ± /N±_ ~ 0.062. This result is in good agreement with 
the predictions of percolation theory. 
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In this paper we show that measuring the distributions of the electric field of TM modes 
of a 3D chaotic rough cavity with the translational symmetry one can find the dependence 
of the number of nodal domains Ktv ± lying on the cross-sectional planes of the cavity on the 
two-dimensional level number iVj_. The translational symmetry means that the cross-section 
of the billiard is invariant under translation along z direction. 

In the seminal papers Blum et al. [l| and Bogomolny and Schmit [2] showed that the 
distributions of the number of nodal domains in two-dimensional (2D) systems can be used to 
distinguish between the systems with integrable and chaotic underlying classical dynamics. 
The theoretical findings have been tested in a series of experiments with chaotic microwave 
2D rough billiards 3], |j, 0] . 

Due to severe experimental problems there are very few experimental studies devoted 
;o 3D chaotic microwave cavities jf], Q, [§], 9], [hJ. In a pioneering experiment Deus et al. 
(| have been measured eigenfrequencies of the 3D chaotic (irregular) microwave cavity 
in order to confirm that their distribution displays behavior characteristic for classically 
chaotic quantum systems, viz., the Wigner distribution. In other important experiments 
the periodic orbits [7], the distributions and the correlation function of the frequency shifts 
caused by the external perturbation [s, 9] and a trace formula for chaotic 3D cavities [h]] 
have been respectively studied. Quite recently the spatial correlation functions of the 3D 
experimental microwave chaotic rough billiard with the translational symmetry have been 



studied by Tymoshchuk et al 



111 ]. Three-dimensional chaotic cavities and properties of 



random electromagnetic vector field have been also studied in several theoretical papers 

The important feature of 3D cavities with the translational symmetry is connected with 
the fact that their modes can be classified into transverse electric (TE) and transverse mag- 
netic (TM). Although, there is no analogy between quantum billiards and electromagnetic 
cavities in three dimensions, the TM modes are especially important because they allow for 
the simulation of 2D quantum billiards on cross-sectional planes of 3D cavities. 

In the experiment we used 3D cavity with the translational symmetry in the shape of 
a rough half-circle (Fig. [1]) with the height h = 60 mm. The cavity was made of polished 
aluminium. The upper and bottom walls of the cavity were attached to the sidewalls with 
48 screws in order to make good electrical contact. 

Assuming that the direction of the translational symmetry of the cavity is along the z- 
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FIG. 1: Sketch of the chaotic half-circular 3D microwave rough billiard in the xy plane. Dimensions 
are given in cm. The cavity sidewalls are marked by 1 and 2 (see text). Squared wave functions 
\^N,p{Rc, 0)\ 2 were evaluated on a half-circle of fixed radius R c = 9.25 cm. Billiard's rough bound- 
ary T is marked with the bold line. The white circle centered at x = 8.12 cm and y = 4.13 cm marks 
the position of the hole drilled in the upper wall of the cavity. The hole was used to introduce the 
perturber inside the cavity in order to measure the z-component of the electric field distributions 



axis the boundary conditions at z = and z = h demand that the z dependence of the 
^-component of the electric and magnetic fields £jv,p(x) and -Bjv, p (x) of TM modes be in 
the form £at, p (x) = E NtP (x,y,z) = A N ^ N , p (x,y)f p (z), where f p (z) = cos(pnz/h), p = 
0, 1, 2 . . ., An, p is the normalization constant and _Bat !P (x) = 0. The dependence of En jP (x) 
on the plane cross section coordinates we denote by the amplitude ipN :P (x,y) = E N;P (x,y). 
Then, the amplitude ipN, P (x,y) satisfies the Helmholtz equation 



where is two-dimensional Laplacian operator and /c/v, p = {k 2 N — (p-K/h) 2 ) l l 2 is the effective 
wave vector. The wave vector = 2711'n/c, where v N is the resonance frequency of the 
level N and c is the speed of light in the vacuum. One can easily see that the equation 
(1) is equivalent to the Schrodinger equation (in units H — 1) describing a particle of mass 



(A ± + k 2 Np )ipN, P ( x ^y) = °> 



(1) 
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m = 1/2 with the kinetic energy k 2 N in an external potential V = (pn/h) 2 [15]. Therefore, 
3D microwave cavities can be effectively used beyond the standard 2D frequency limit (the 



case p = 0) [16(| in simulation of quantum systems. The amplitude tp^^x, y) fulfills Dirichlet 
boundary conditions on the sidewalls of the billiard and therefore, throughout the text it is 
also called the wave functions iPn{%, v)- It is worth noting that the full electric field E^ jP (x.) 
satisfies Neumann boundary conditions at the top and the bottom of the cavity. 

The measurements of £^at jP (x) of a 3D microwave cavity allowed us to test experimentally 
an important finding of the papers by Blum et al. [l| and Bogomolny and Schmit which 
connects the number of nodal domains of 2D billiards with the level number N. We will 
show that for the 3D cavities with the translational symmetry the number of nodal domains 
Kjy ± lying on the cross-sectional planes of the cavity is connected with the two-dimensional 
level number N±. The condition E^^{-x)\ z=const = on the cross-sectional planes of the 
cavity determines a set of nodal lines which separate regions (nodal domains) with opposite 
signs of the electric field distribution -E7v iP (x)| 2=cons t. 

The value of the level number N of the 3D cavity was evaluated from the Balian-Bloch 
formula [171 ] . 

where k is the wave vector, V = (9.43 ± 0.01) ■ 10~ 4 m 3 is the volume of the cavity and 
J s ^ = 0.932 m ±0.005 m is the surface curvature averaged over the surface of the cavity. 
We used this formula because of the relatively low quality factor of the cavity (Q ~ 4000) 
some resonances overlapped. 

The two-dimensional level number N±_ is defined by the standard Weyl-Bloch formula 
N ± = ^k 2 Np - £k NtP , where A = (1.572 ± 0.002) • 10~ 2 m 2 and P = 0.537 m ±0.001 m are 
the cross-sectional plane area of the cavity and its perimeter, respectively. 

The cavity sidewalls consist of two segments (see Fig.[T]). The rough segment 1 is described 
on the cross-sectional planes by the radius function R(6) = R + Dm=2 a m sm(m9 + 4> m ), 
where the mean radius i?o=10.0 cm, M = 20, a m and (f) m are uniformly distributed on 
[0.084,0.091] cm and [0,27r], respectively, and < 6 < it. (For convenience, the polar 
coordinates r and 9 are used instead of the Cartesian ones x and y.) It is worth noting 
that following our earlier experience we decided to use a rough desymmetrized half- 

circular cavity instead of a rough circular cavity, because the first one lowers the number 
of nearly degenerated eigenfrequencies. Additionally, a half-circular geometry of the cavity 



was suitable in the procedure of accurate measurements of the electric field distributions 
inside the billiard. 

The roughness of a billiard on the cross-sectional planes can be characterized by the 
function k(9) = (dR/d6)/ Rq. For our microwave billiard we have the angle average k = 
_ik 2 (6)} e ) 1/2 ~ 0.400. The value of k is much above the chaos border k c = M" 5 / 2 = 0.00056 



201 ] which indicates that in such a billiard the classical dynamics is diffusive in orbital 



momentum due to collisions with the r oug h boundary. 

The other properties of the billiard 2l| are also determined by the roughness parameter 
k. The amplitudes i/jN,p{r, 9) are localized for the two-dimensional level number N± < N e = 
1/128/c 4 . Because of a large value of the roughness parameter k the localization border lies 
very low, N e ~ 1. The border of Breit-Wigner regime is N\y = M 2 /48k 2 ~ 52. It means 
that between N e < N± < N\y Wigner ergodicity [21J ought to be observed and for N± > Nw 
Shnirelman ergodicity should emerge. 

In order to measure the amplitudes if)N&{ r , 9) of the 3D electric field distributions we used 
an effective method described in 22]. It is based on the perturbation technique and prepa- 
ration of the "trial functions". In this method the amplitudes ipN^, 9) (electric field distri- 
bution Epj(r, 9) inside the cavity) are determined from the form of electric field E^^ p (R c , 9) 
evaluated on a half-circle of fixed radius R c (see Fig. [1]). The first step in evaluation of 



En, p (R c ,9) is measurement of \ E N J R r _, 9)\ 2 . The perturbation technique developed in [23] 



and used successfully in 23|, |24j, |25|, [26] was implemented for this purpose. In this method 
a small perturber is introduced inside the cavity to alter its resonant frequency. 

The perturber (4.0 mm in length and 0.3 mm in diameter, oriented in z-direction) was 
moved by the stepper motor via the Kevlar line hidden in the groove (0.4 mm wide, 1.0 mm 
deep) made in the cavity's bottom wall along the half-circle R c . Before closing the cavity 
we carefully inspected whether the pin moves smoothly, oriented in vertical position. Using 
such a perturber we had no positive frequency shifts that would exceed the uncertainty of 
frequency shift measurements (15 kHz). 

In order to determine the dependence of the electric field distributions En >p (x) on the z 
coordinate and to estimate the wave vector fc 3 = pw/h we measured the electric field inside 
the 3D cavity along the z-axis. The perturber (4.5 mm in length and 0.3 mm in diameter) 
was attached to the Kevlar line and moved by the stepper motor. It entered and exited the 
cavity by small holes (0.4 mm) drilled in the upper and the bottom walls of the cavity. The 
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FIG. 2: Panel (a): Squared wave function 1^430,2(^0 &)\ 2 (in arbitrary units) measured on a half- 
circle with radius R c = 9.25 cm (^430 — 11.50 GHz). Panel (b): Squared z-component of the 
electric field distribution |/2(z)| 2 measured at r = 9.11 cm and 9 = 0.47 radians. 



both holes were located at the position: r = 9.11 cm, 6 = 0.47 radians. 

To eliminate the variation of resonant frequencies connected with the thermal expansion 
of the aluminium cavity the temperature of the cavity was stabilized with the accuracy of 
0.05 deg. 

Using a field perturbation technique we were able to measure squared wave functions 
\ip N ,p( R ci °)\ 2 for 80 ™ modes within the region 2 < N < 489. The range of corresponding 
eigenfrequencies was from u 2 ~ 2.47 GHz to z/ 489 ~ 11.99 GHz. The measurements were 
performed at 0.36 mm steps along a half-circle with fixed radius R c = 9.25 cm. This step 
was small enough to reveal in details the space structure of high-lying levels. 
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In Fig. [2](a) and Fig. 121(b) we show the examples of the squared amplitude \ipN, P {Rc, #)| 2 
and the squared ^-component of the field, respectively, evaluated for the level number N = 
430. 

The perturbation method used in our measurements allows us to extract information 
about the modulus of the wave function amplitude \ipN,p(R c ,@)\ at any given point of the 
cross-sectional plane z = but it doesn't allow to determine the sign of ipN, P {R c , 0). In order 
to obtain information about the sign oi iI)m :P (R c , 6) we used the method of the "trial wave 



function" precisely described in 



St, 



221. 



The amplitudes ipN,p(j", 6) of the electric field distributions of a rough half-circular 3D 
billiard may be expanded in terms of circular waves (here only odd states in expansion are 
considered) 

L 

ipN, P {r, 0) = a sJ s (kN, P r) sin(s0), (5) 

8=1 

where J s is the Bessel function of order s. 

In Eq. (5) the number of basis functions is limited to L — k N ^ p r max = l^ ax , where 
Tmax = 10.64 cm is the maximum radius of the cavity. l^ ax = k^^ v r max is a semiclassical 
estimate for the maximum possible angular momentum for a given k^. Circular waves with 
angular momentum s > L correspond to evanescent waves and can be neglected. Coefficients 
a s may be extracted from the "trial wave function" i/jn,p(Rc, 0) via 



Js{k N , P Rc)V f *Pn, p ( R ci 0) sm(s9)de. (6) 

J 



Due to experimental uncertainties and the finite step size in the measurements of 
\ipN,p(Rc, 9)\ 2 the wave functions ipN, P (r,9) are not exactly zero at the boundary T. As 
the quantitative measure of the sign assignment quality we chose the integral / = 
7 It l'0JV,p( r j 9)\ 2 dl calculated along the billiard's rough boundary T, where 7 is length of 
T. For correctly reconstructed wave functions the integral I was several times smaller than 
in the case of not correctly reconstructed ones. 

It is worth noting that since the pin is attached to the line it cannot be stuck. However, 
one may assume that during the movement the pin may be accidentally, from time to 
time, slightly slanted, adding small " a noise-like component" to the measured electric field. 
The formula (5) shows that each wave function is expanded in terms of L circular waves 
which filters out noise-like higher frequency Fourier components from the reconstructed wave 
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FIG. 3: The "trial wave function" ^430 i(R C i Q) (hi arbitrary units) with the correctly assigned 
signs, which was used in the reconstruction of the wave function £430,2 ( r > 0, z) of the billiard (see 

Fig. m). 

function. The same filtering removes out the influence of the experimental uncertainties of 
frequency shifts on the reconstructed wave functions. 

In Fig. [3] we show the "trial wave function" ^430,2 {Rc, 9) with the correctly assigned signs, 
which was used in the reconstruction of the wave function ^430,2 (f) 0) of the billiard (see Fig. 
HJ) . In Fig. H] different nodal domains are separated by the bold full lines. 

Using the method of the "trial wave function" we were able to reconstruct 75 experi- 
mental wave functions ^jv jP (r, 0), which belonged to TM modes of the rough half-circular 
3D billiard with the level number N between 2 and 489. The remaining wave functions 
belonging to TM modes, from the range N = 2 — 489, were not reconstructed because of 
near-degeneration of the neighboring eigenfrequencies or due to the problems with the mea- 
surements of \ipN tP (R c ,8)\ 2 along a half-circle coinciding for its significant part with one of 
the nodal lines of ipN, P ( r , 0). 

The borders of Breit-Wigner and Shnirelman ergodicities are not sharp. Therefore, to 
check ergodicity of the billiard's wave functions i/;N jP (r, @), especially close to the borders, one 
should use some additional measures such as e.g., calculation of the structures of their energy 
surfaces [20]. For this reason we extracted wave function amplitudes = (n,l\N,p) in 
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FIG. 4: The reconstructed wave function ^430,2 ( r ,6) of the chaotic half-circular microwave rough 
billiard. The amplitudes have been converted into a grey scale with white corresponding to large 
positive and black corresponding to large negative values, respectively. The structure of the nodal 
lines are shown by the bold full lines. Dimensions of the billiard are given in cm. In the figure 
the z dependence of the electric field distribution ^430,2^, 0, z)\x=o cx ^430,2 ('"j ^) |x=o/2 (^) is also 
shown. 

the basis n, I of a half-circular billiard with radius r max , where n — 1, 2, 3 . . . enumerates the 
zeros of the Bessel functions and I — 1, 2, 3 ... is the angular quantum number. The moduli 
of amplitudes \C^' P '\ and their projections into the energy surface for the experimental 
wave functions ?/>46o,o( r ; Q), ^430,2 ( r , 0) and ^463,4 ( r ; 9) are shown in Fig. [5la-c). As expected, 
on the border of the regimes of Breit-Wigner and Shnirelman ergodicity the wave functions 
■0460,0 ( r , 6) {N± = 65) and ^430,2 {f, 0) (N± = 50) are extended homogeneously over the 
whole energy surface The wave function ^463,4 ( r > $)> -^Y± — 16, which lies closer to 

the localization boarder, is also extended along the energy surface, however it displays the 
tendency to localization in n,l basis (see Fig. E](c)). The full lines on the projection planes 
in Fig. [5](a-c) mark the energy surface of a half-circular billiard H(n,l) = k 2 Np estimated 
from the semiclassical formula [19[: J {l^ ax ) 2 — I 2 — I arctan(/ _1 y / {l 1 ^^) 2 — I 2 ) + 7r/4 = irn. 
It is clearly visible that the peaks | are spread almost perfectly along the lines marking 

the energy surface. 

The number of nodal domains Kjv ± on the cross-sectional plane z = vs. the level 
number N± in the chaotic 3D microwave rough billiard is plotted in Fig. [61 The full line in 
Fig. [6] shows the least squares fit Kjv x = a\N± + biy/N± of the experimental data, where 
ai = 0.059 ±0.029, 61 = 0.991 ±0.190. The coefficient a x = 0.059 ±0.029 coincides with the 
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FIG. 5: Structure of the energy surface of the wave functions lying close to the boarder of the 
regimes of Breit-Wigner and Shnirelman ergodicity (Nw = 52), panels (a) and (b), and for the 
low wave function in the regime of Breit-Wigner ergodicity, panel (c). Panel (a): The moduli of 
amplitudes IC^ 60 ' - 1 ! for the wave function ^46o,o( r > 0), N± = 65, lying in the regime of Shnirelman 
ergodicity. Panel (b): The moduli of amplitudes |C^ 30,2 ' ) | for the wave function ^430,2 (?", &), N± = 
50, in the regime of Breit-Wigner ergodicity . Panel (c): The moduli of amplitudes |C^t 63 ' 4 ^| for the 
wave function ^463,4 (?*> lying in the regime of Breit-Wigner ergodicity close to the localization 
boarder. Full lines show the semiclassical estinlDtion of the energy surface (see text). 
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FIG. 6: The number of nodal domains Nat ± (full circles) on the cross-section planes of the chaotic 
half-circular 3D microwave rough billiard. Full line shows the least squares fit ^at ± = a\N±+biy/ N± 
to the experimental data (see text), where a% = 0.059 ± 0.029, b\ = 0.991 ± 0.190. The prediction 
of the theory of Bogomolny and Schmit [2] a\ = 0.062. 

prediction of the percolation model of Bogomolny and Schmit j^j Hn ± /N± ~ 0.062 within 
the error limits. The relatively large uncertainty of the coefficient cq is connected with the 
fact that in the least squares fit procedure we used only 27 higher states with N± > 20. The 
states with lower N± were not taken into account because they were not fully chaotic (see Fig. 
IHc)). The second term in the least squares fit corresponds to a contribution of boundary 
domains, i.e. domains, which include the billiard boundary. Numerical calculations of Blum 
et al. [l| performed for the Sinai and stadium billiards showed that the number of boundary 
domains scales as the number of the boundary intersections, that is as y/N±. Our results 
clearly suggest that in the rough billiard, at the level numbers 20 < N± < 65, the boundary 
domains also significantly influence the scaling of the number of nodal domains Kat ± , leading 
to the departure from the predicted scaling Kat ± ~ N±. 

In summary, we measured the wave functions of the chaotic 3D rough microwave billiard 
with the translational symmetry up to the level number N = 489. We showed that for the 
two-dimensional level numbers 20 < iV^ < 65 the scaling of the number of nodal domains 
Kjy ± significantly departures from the predicted scaling #, N± ~ N±, which suggests that the 



11 



boundary domains influence the scaling [2|. In the limit N± — > oo the least squares fit of 
the experimental data yields the asymptotic number of nodal domains ^n ± /N± ~ a\ — 
0.059 ±0.029 that is close to the theoretical prediction ^n ± /N± ~ 0.062. Finally, our results 
show that 3D microwave cavities with the translational symmetry can be effectively used 
beyond the standard 2D frequency limit in simulation of quantum systems. 
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